Exploration of the potential common pathogenic mechanisms in COVID-19 and silicosis by using bioinformatics and system biology.

Silicosis is an occupational lung disease that is common worldwide. In recent years, coronavirus disease 2019 (COVID-19) has provided daunting challenges to public healthcare systems globally. Although multiple studies have shown a close link between COVID-19 and other respiratory diseases, the inter-relational mechanisms between COVID-19 and silicosis remain unclear. This study aimed to explore the shared molecular mechanisms and drug targets of COVID-19 and silicosis. Gene expression profiling identified four modules that were most closely associated with both diseases. Furthermore, we performed functional analysis and constructed a protein-protein interaction network. Seven hub genes (budding uninhibited by benzimidazoles 1 [BUB1], protein regulator of cytokinesis 1 [PRC1], kinesin family member C1 [KIFC1], ribonucleotide reductase regulatory subunit M2 [RRM2], cyclin-dependent kinase inhibitor 3 [CDKN3], Cyclin B2 [CCNB2], and minichromosome maintenance complex component 6 [MCM6]) were involved in the interaction between COVID-19 and silicosis. We investigated how diverse microRNAs and transcription factors regulate these seven genes. Subsequently, the correlation between the hub genes and infiltrating immune cells was explored. Further in-depth analyses were performed based on single-cell transcriptomic data from COVID-19, and the expression of hub-shared genes was characterized and located in multiple cell clusters. Finally, molecular docking results reveal small molecular compounds that may improve COVID-19 and silicosis. The current study reveals the common pathogenesis of COVID-19 and silicosis, which may provide a novel reference for further research.

SEMtree: tree-based structure learning methods with structural equation models.

MOTIVATION: With the exponential growth of expression and protein-protein interaction (PPI) data, the identification of functional modules in PPI networks that show striking changes in molecular activity or phenotypic signatures becomes of particular interest to reveal process-specific information that is correlated with cellular or disease states. This requires both the identification of network nodes with reliability scores and the availability of an efficient technique to locate the network regions with the highest scores. In the literature, a number of heuristic methods have been suggested. We propose SEMtree(), a set of tree-based structure discovery algorithms, combining graph and statistically interpretable parameters together with a user-friendly R package based on structural equation models framework. RESULTS: Condition-specific changes from differential expression and gene-gene co-expression are recovered with statistical testing of node, directed edge, and directed path difference between groups. In the end, from a list of seed (i.e. disease) genes or gene P-values, the perturbed modules with undirected edges are generated with five state-of-the-art active subnetwork detection methods. The latter are supplied to causal additive trees based on Chu-Liu-Edmonds' algorithm (Chow and Liu, Approximating discrete probability distributions with dependence trees. IEEE Trans Inform Theory 1968;14:462-7) in SEMtree() to be converted in directed trees. This conversion allows to compare the methods in terms of directed active subnetworks. We applied SEMtree() to both Coronavirus disease (COVID-19) RNA-seq dataset (GEO accession: GSE172114) and simulated datasets with various differential expression patterns. Compared to existing methods, SEMtree() is able to capture biologically relevant subnetworks with simple visualization of directed paths, good perturbation extraction, and classifier performance. AVAILABILITY AND IMPLEMENTATION: SEMtree() function is implemented in the R package SEMgraph, easily available at https://CRAN.R-project.org/package=SEMgraph.

A robust protein language model for SARS-CoV-2 protein-protein interaction network prediction.

Protein-protein interaction is one of the ways viruses interact with their hosts. Therefore, identifying protein interactions between viruses and hosts helps explain how virus proteins work, how they replicate, and how they cause disease. SARS-CoV-2 is a new type of virus that emerged from the coronavirus family in 2019 and caused a worldwide pandemic. Detection of human proteins interacting with this novel virus strain plays an important role in monitoring the cellular process of virus-associated infection. Within the scope of the study, a natural language processing-based collective learning method is proposed for the prediction of potential SARS-CoV-2-human PPIs. Protein language models were obtained with the prediction-based word2Vec and doc2Vec embedding methods and the frequency-based tf-idf method. Known interactions were represented by proposed language models and traditional feature extraction methods (conjoint triad and repeat pattern), and their performances were compared. The interaction data were trained with support vector machine, artificial neural network (ANN), k-nearest neighbor (KNN), naive Bayes (NB), decision tree (DT), and ensemble algorithms. Experimental results show that protein language models are a promising protein representation method for protein-protein interaction prediction. The term frequency-inverse document frequency-based language model performed the SARS-CoV-2 protein-protein interaction estimation with an error of 1.4%. Additionally, the decisions of high-performing learning models for different feature extraction methods were combined with a collective voting approach to make new interaction predictions. For 10,000 human proteins, 285 new potential interactions were predicted, with models combining decisions.

De novo design of protein interactions with learned surface fingerprints.

Physical interactions between proteins are essential for most biological processes governing life1. However, the molecular determinants of such interactions have been challenging to understand, even as genomic, proteomic and structural data increase. This knowledge gap has been a major obstacle for the comprehensive understanding of cellular protein-protein interaction networks and for the de novo design of protein binders that are crucial for synthetic biology and translational applications2-9. Here we use a geometric deep-learning framework operating on protein surfaces that generates fingerprints to describe geometric and chemical features that are critical to drive protein-protein interactions10. We hypothesized that these fingerprints capture the key aspects of molecular recognition that represent a new paradigm in the computational design of novel protein interactions. As a proof of principle, we computationally designed several de novo protein binders to engage four protein targets: SARS-CoV-2 spike, PD-1, PD-L1 and CTLA-4. Several designs were experimentally optimized, whereas others were generated purely in silico, reaching nanomolar affinity with structural and mutational characterization showing highly accurate predictions. Overall, our surface-centric approach captures the physical and chemical determinants of molecular recognition, enabling an approach for the de novo design of protein interactions and, more broadly, of artificial proteins with function.

Bioinformatics Analysis of Candidate Genes and Pathways Related to Hepatocellular Carcinoma in China: A Study Based on Public Databases.

Background and Objective: Hepatocellular carcinoma (HCC) is a highly aggressive malignant tumor of the digestive system worldwide. Chronic hepatitis B virus (HBV) infection and aflatoxin exposure are predominant causes of HCC in China, whereas hepatitis C virus (HCV) infection and alcohol intake are likely the main risk factors in other countries. It is an unmet need to recognize the underlying molecular mechanisms of HCC in China. Methods: In this study, microarray datasets (GSE84005, GSE84402, GSE101685, and GSE115018) derived from Gene Expression Omnibus (GEO) database were analyzed to obtain the common differentially expressed genes (DEGs) by R software. Moreover, the gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed by using Database for Annotation, Visualization and Integrated Discovery (DAVID). Furthermore, the protein-protein interaction (PPI) network was constructed, and hub genes were identified by the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape, respectively. The hub genes were verified using Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and Kaplan-Meier Plotter online databases were performed on the TCGA HCC dataset. Moreover, the Human Protein Atlas (HPA) database was used to verify candidate genes' protein expression levels. Results: A total of 293 common DEGs were screened, including 103 up-regulated genes and 190 down-regulated genes. Moreover, GO analysis implied that common DEGs were mainly involved in the oxidation-reduction process, cytosol, and protein binding. KEGG pathway enrichment analysis presented that common DEGs were mainly enriched in metabolic pathways, complement and coagulation cascades, cell cycle, p53 signaling pathway, and tryptophan metabolism. In the PPI network, three subnetworks with high scores were detected using the Molecular Complex Detection (MCODE) plugin. The top 10 hub genes identified were CDK1, CCNB1, AURKA, CCNA2, KIF11, BUB1B, TOP2A, TPX2, HMMR and CDC45. The other public databases confirmed that high expression of the aforementioned genes related to poor overall survival among patients with HCC. Conclusion: This study primarily identified candidate genes and pathways involved in the underlying mechanisms of Chinese HCC, which is supposed to provide new targets for the diagnosis and treatment of HCC in China.

Functional Analysis of Monkeypox and Interrelationship between Monkeypox and COVID-19 by Bioinformatic Analysis.

The outbreak of monkeypox may be considered a novel and urgent threat after the coronavirus disease (COVID-19). No wide-ranging studies have been conducted on this disease since it was first reported. We systematically assessed the functional role of gene expression in cells infected with the monkeypox virus using transcriptome profiling and compared the functional relation with that of COVID-19. Based on the Gene Expression Omnibus database, we obtained 212 differentially expressed genes (DEGs) of GSE36854 and GSE21001 of monkeypox datasets. Enrichment analyses, including KEGG and gene ontology (GO) analyses, were performed to identify the common function of 212 DEGs of GSE36854 and GSE21001. CytoHubba and Molecular Complex Detection were performed to determine the core genes after a protein-protein interaction (PPI). Metascape/COVID-19 was used to compare DEGs of monkeypox and COVID-19. GO analysis of 212 DEGs of GSE36854 and GSE21001 for monkeypox infection showed cellular response to cytokine stimulus, cell activation, and cell differentiation regulation. KEGG analysis of 212 DEGs of GSE36854 and GSE21001 for monkeypox infection showed involvement of monkeypox in COVID-19, cytokine-cytokine receptor interaction, inflammatory bowel disease, atherosclerosis, TNF signaling, and T cell receptor signaling. By comparing our data with published transcriptome of severe acute respiratory syndrome coronavirus 2 infections in other cell lines, the common function of monkeypox and COVID-19 includes cytokine signaling in the immune system, TNF signaling, and MAPK cascade regulation. Thus, our data suggest that the molecular connections identified between COVID-19 and monkeypox elucidate the causes of monkeypox.

Exploration of the Shared Molecular Mechanisms between COVID-19 and Neurodegenerative Diseases through Bioinformatic Analysis.

The COVID-19 pandemic has caused millions of deaths and remains a major public health burden worldwide. Previous studies found that a large number of COVID-19 patients and survivors developed neurological symptoms and might be at high risk of neurodegenerative diseases, such as Alzheimer's disease (AD) and Parkinson's disease (PD). We aimed to explore the shared pathways between COVID-19, AD, and PD by using bioinformatic analysis to reveal potential mechanisms, which may explain the neurological symptoms and degeneration of brain that occur in COVID-19 patients, and to provide early intervention. In this study, gene expression datasets of the frontal cortex were employed to detect common differentially expressed genes (DEGs) of COVID-19, AD, and PD. A total of 52 common DEGs were then examined using functional annotation, protein-protein interaction (PPI) construction, candidate drug identification, and regulatory network analysis. We found that the involvement of the synaptic vesicle cycle and down-regulation of synapses were shared by these three diseases, suggesting that synaptic dysfunction might contribute to the onset and progress of neurodegenerative diseases caused by COVID-19. Five hub genes and one key module were obtained from the PPI network. Moreover, 5 drugs and 42 transcription factors (TFs) were also identified on the datasets. In conclusion, the results of our study provide new insights and directions for follow-up studies of the relationship between COVID-19 and neurodegenerative diseases. The hub genes and potential drugs we identified may provide promising treatment strategies to prevent COVID-19 patients from developing these disorders.

Computational network analysis of host genetic risk variants of severe COVID-19.

BACKGROUND: Genome-wide association studies have identified numerous human host genetic risk variants that play a substantial role in the host immune response to SARS-CoV-2. Although these genetic risk variants significantly increase the severity of COVID-19, their influence on body systems is poorly understood. Therefore, we aim to interpret the biological mechanisms and pathways associated with the genetic risk factors and immune responses in severe COVID-19. We perform a deep analysis of previously identified risk variants and infer the hidden interactions between their molecular networks through disease mapping and the similarity of the molecular functions between constructed networks. RESULTS: We designed a four-stage computational workflow for systematic genetic analysis of the risk variants. We integrated the molecular profiles of the risk factors with associated diseases, then constructed protein-protein interaction networks. We identified 24 protein-protein interaction networks with 939 interactions derived from 109 filtered risk variants in 60 risk genes and 56 proteins. The majority of molecular functions, interactions and pathways are involved in immune responses; several interactions and pathways are related to the metabolic and cardiovascular systems, which could lead to multi-organ complications and dysfunction. CONCLUSIONS: This study highlights the importance of analyzing molecular interactions and pathways to understand the heterogeneous susceptibility of the host immune response to SARS-CoV-2. We propose new insights into pathogenicity analysis of infections by including genetic risk information as essential factors to predict future complications during and after infection. This approach may assist more precise clinical decisions and accurate treatment plans to reduce COVID-19 complications.

Identification of differentially expressed genes and pathways in kidney of ANCA-associated vasculitis by integrated bioinformatics analysis.

The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) is a systematic of relatively rare autoimmune diseases with unknown cause. Kidney involvement is one of the most common clinical manifestations, and the degree of renal damage is closely associated with the development and prognosis of AAV. In this study, we utilized the Robust Rank Aggreg (RRA) method in R to integrate GSE104948, GSE104954, GSE108109, GSE108112, and GSE108113 profile datasets loaded from Gene Expression Omnibus (GEO) database and identified a set of differentially expressed genes (DEGs) in kidney between AAV patients and living donors. Then, the results of gene ontology (GO) functional annotation showed that immunity and metabolism involved process of AAV both in glomerulus and tubulointerstitial. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that following pathways, such as complement and coagulation cascades pathway; Staphylococcus aureus infection; disease-COVID-19; and systemic lupus erythematosus (SLE) pathway play a crucial role in AAV. Next, the results analyzed by protein-protein interaction (PPI) network and Cytoscape software exhibited the hub genes ALB, TYROBP, and CYBB existed in both glomerular and tubulointerstitial compartments datasets. Finally, KEGG analysis using genes of two most important modules also further validated complement and coagulation cascades pathway and S. aureus infection existed both in glomerulus and tubulointerstitial compartments datasets. In conclusion, this study identified key genes and pathways involved in kidney of AAV, which was benefit to further uncover the mechanisms underlying the development and progress of AAV, biomarkers, and potential therapeutic targets as well.

Peripheral Blood Genes Crosstalk between COVID-19 and Sepsis.

Severe coronavirus disease 2019 (COVID-19) has led to a rapid increase in death rates all over the world. Sepsis is a life-threatening disease associated with a dysregulated host immune response. It has been shown that COVID-19 shares many similarities with sepsis in many aspects. However, the molecular mechanisms underlying sepsis and COVID-19 are not well understood. The aim of this study was to identify common transcriptional signatures, regulators, and pathways between COVID-19 and sepsis, which may provide a new direction for the treatment of COVID-19 and sepsis. First, COVID-19 blood gene expression profile (GSE179850) data and sepsis blood expression profile (GSE134347) data were obtained from GEO. Then, we intersected the differentially expressed genes (DEG) from these two datasets to obtain common DEGs. Finally, the common DEGs were used for functional enrichment analysis, transcription factor and miRNA prediction, pathway analysis, and candidate drug analysis. A total of 307 common DEGs were identified between the sepsis and COVID-19 datasets. Protein-protein interactions (PPIs) were constructed using the STRING database. Subsequently, hub genes were identified based on PPI networks. In addition, we performed GO functional analysis and KEGG pathway analysis of common DEGs, and found a common association between sepsis and COVID-19. Finally, we identified transcription factor-gene interaction, DEGs-miRNA co-regulatory networks, and protein-drug interaction, respectively. Through ROC analysis, we identified 10 central hub genes as potential biomarkers. In this study, we identified SARS-CoV-2 infection as a high risk factor for sepsis. Our study may provide a potential therapeutic direction for the treatment of COVID-19 patients suffering from sepsis.

Recognition of Differentially Expressed Molecular Signatures and Pathways Associated with COVID-19 Poor Prognosis in Glioblastoma Patients.

Glioblastoma (GBM) is a type of brain cancer that is typically very aggressive and difficult to treat. Glioblastoma cases have been reported to have increased during COVID-19. The mechanisms underlying this comorbidity, including genomic interactions, tumor differentiation, immune responses, and host defense, are not completely explained. Therefore, we intended to investigate the differentially expressed shared genes and therapeutic agents which are significant for these conditions by using in silico approaches. Gene expression datasets of GSE68848, GSE169158, and GSE4290 studies were collected and analyzed to identify the DEGs between the diseased and the control samples. Then, the ontology of the genes and the metabolic pathway enrichment analysis were carried out for the classified samples based on expression values. Protein-protein interactions (PPI) map were performed by STRING and fine-tuned by Cytoscape to screen the enriched gene module. In addition, the connectivity map was used for the prediction of potential drugs. As a result, 154 overexpressed and 234 under-expressed genes were identified as common DEGs. These genes were found to be significantly enriched in the pathways involved in viral diseases, NOD-like receptor signaling pathway, the cGMP-PKG signaling pathway, growth hormone synthesis, secretion, and action, the immune system, interferon signaling, and the neuronal system. STAT1, CXCL10, and SAMDL were screened out as the top 03 out of the top 10 most critical genes among the DEGs from the PPI network. AZD-8055, methotrexate, and ruxolitinib were predicted to be the possible agents for the treatment. The current study identified significant key genes, common metabolic signaling networks, and therapeutic agents to improve our perception of the common mechanisms of GBM-COVID-19.

A network view of human immune system and virus-human interaction.

The immune system is highly networked and complex, which is continuously changing as encountering old and new pathogens. However, reductionism-based researches do not give a systematic understanding of the molecular mechanism of the immune response and viral pathogenesis. Here, we present HUMPPI-2022, a high-quality human protein-protein interaction (PPI) network, containing > 11,000 protein-coding genes with > 78,000 interactions. The network topology and functional characteristics analyses of the immune-related genes (IRGs) reveal that IRGs are mostly located in the center of the network and link genes of diverse biological processes, which may reflect the gene pleiotropy phenomenon. Moreover, the virus-human interactions reveal that pan-viral targets are mostly hubs, located in the center of the network and enriched in fundamental biological processes, but not for coronavirus. Finally, gene age effect was analyzed from the view of the host network for IRGs and virally-targeted genes (VTGs) during evolution, with IRGs gradually became hubs and integrated into host network through bridging functionally differentiated modules. Briefly, HUMPPI-2022 serves as a valuable resource for gaining a better understanding of the composition and evolution of human immune system, as well as the pathogenesis of viruses.

Discovery of host-directed modulators of virus infection by probing the SARS-CoV-2-host protein-protein interaction network.

The ongoing coronavirus disease 2019 (COVID-19) pandemic has highlighted the need to better understand virus-host interactions. We developed a network-based method that expands the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-host protein interaction network and identifies host targets that modulate viral infection. To disrupt the SARS-CoV-2 interactome, we systematically probed for potent compounds that selectively target the identified host proteins with high expression in cells relevant to COVID-19. We experimentally tested seven chemical inhibitors of the identified host proteins for modulation of SARS-CoV-2 infection in human cells that express ACE2 and TMPRSS2. Inhibition of the epigenetic regulators bromodomain-containing protein 4 (BRD4) and histone deacetylase 2 (HDAC2), along with ubiquitin-specific peptidase (USP10), enhanced SARS-CoV-2 infection. Such proviral effect was observed upon treatment with compounds JQ1, vorinostat, romidepsin and spautin-1, when measured by cytopathic effect and validated by viral RNA assays, suggesting that the host proteins HDAC2, BRD4 and USP10 have antiviral functions. We observed marked differences in antiviral effects across cell lines, which may have consequences for identification of selective modulators of viral infection or potential antiviral therapeutics. While network-based approaches enable systematic identification of host targets and selective compounds that may modulate the SARS-CoV-2 interactome, further developments are warranted to increase their accuracy and cell-context specificity.

A protein-protein interaction map reveals that the Coxiella burnetii effector CirB inhibits host proteasome activity.

Coxiella burnetii is the etiological agent of the zoonotic disease Q fever, which is featured by its ability to replicate in acid vacuoles resembling the lysosomal network. One key virulence determinant of C. burnetii is the Dot/Icm system that transfers more than 150 effector proteins into host cells. These effectors function to construct the lysosome-like compartment permissive for bacterial replication, but the functions of most of these effectors remain elusive. In this study, we used an affinity tag purification mass spectrometry (AP-MS) approach to generate a C. burnetii-human protein-protein interaction (PPI) map involving 53 C. burnetii effectors and 3480 host proteins. This PPI map revealed that the C. burnetii effector CBU0425 (designated CirB) interacts with most subunits of the 20S core proteasome. We found that ectopically expressed CirB inhibits hydrolytic activity of the proteasome. In addition, overexpression of CirB in C. burnetii caused dramatic inhibition of proteasome activity in host cells, while knocking down CirB expression alleviated such inhibitory effects. Moreover, we showed that a region of CirB that spans residues 91-120 binds to the proteasome subunit PSMB5 (beta 5). Finally, PSMB5 knockdown promotes C. burnetii virulence, highlighting the importance of proteasome activity modulation during the course of C. burnetii infection.

Identification of significant alteration genes, pathways and TFs induced by LPS in ARDS via bioinformatical analysis.

BACKGROUND AND AIMS: Acute respiratory distress syndrome (ARDS) or acute lung injury (ALI) is one of the most common acute thoracopathy with complicated pathogenesis in ICU. The study is to explore the differentially expressed genes (DEGs) in the lung tissue and underlying altering mechanisms in ARDS. METHODS: Gene expression profiles of GSE2411 and GSE130936 were available from GEO database, both of them included in GPL339. Then, an integrated analysis of these genes was performed, including gene ontology (GO) and KEGG pathway enrichment analysis in DAVID database, protein-protein interaction (PPI) network construction evaluated by the online database STRING, Transcription Factors (TFs) forecasting based on the Cytoscape plugin iRegulon, and their expression in varied organs in The Human Protein Atlas. RESULTS: A total of 39 differential expressed genes were screened from the two datasets, including 39 up-regulated genes and 0 down-regulated genes. The up-regulated genes were mainly enriched in the biological process, such as immune system process, innate immune response, inflammatory response, and also involved in some signal pathways, including cytokine-cytokine receptor interaction, Salmonella infection, Legionellosis, Chemokine, and Toll-like receptor signal pathway with an integrated analysis. GBP2, IFIT2 and IFIT3 were identified as hub genes in the lung by PPI network analysis with MCODE plug-in, as well as GO and KEGG re-enrichment. All of the three hub genes were regulated by the predictive common TFs, including STAT1, E2F1, IRF1, IRF2, and IRF9. CONCLUSIONS: This study implied that hub gene GBP2, IFIT2 and IFIT3, which might be regulated by STAT1, E2F1, IRF1, IRF2, or IRF9, played significant roles in ARDS. They could be potential diagnostic or therapeutic targets for ARDS patients.

Comparing protein-protein interaction networks of SARS-CoV-2 and (H1N1) influenza using topological features.

SARS-CoV-2 pandemic first emerged in late 2019 in China. It has since infected more than 298 million individuals and caused over 5 million deaths globally. The identification of essential proteins in a protein-protein interaction network (PPIN) is not only crucial in understanding the process of cellular life but also useful in drug discovery. There are many centrality measures to detect influential nodes in complex networks. Since SARS-CoV-2 and (H1N1) influenza PPINs pose 553 common human proteins. Analyzing influential proteins and comparing these networks together can be an effective step in helping biologists for drug-target prediction. We used 21 centrality measures on SARS-CoV-2 and (H1N1) influenza PPINs to identify essential proteins. We applied principal component analysis and unsupervised machine learning methods to reveal the most informative measures. Appealingly, some measures had a high level of contribution in comparison to others in both PPINs, namely Decay, Residual closeness, Markov, Degree, closeness (Latora), Barycenter, Closeness (Freeman), and Lin centralities. We also investigated some graph theory-based properties like the power law, exponential distribution, and robustness. Both PPINs tended to properties of scale-free networks that expose their nature of heterogeneity. Dimensionality reduction and unsupervised learning methods were so effective to uncover appropriate centrality measures.

Identification of critical genes and molecular pathways in COVID-19 myocarditis and constructing gene regulatory networks by bioinformatic analysis.

BACKGROUND: There is growing evidence of a strong relationship between COVID-19 and myocarditis. However, there are few bioinformatics-based analyses of critical genes and the mechanisms related to COVID-19 Myocarditis. This study aimed to identify critical genes related to COVID-19 Myocarditis by bioinformatic methods, explore the biological mechanisms and gene regulatory networks, and probe related drugs. METHODS: The gene expression data of GSE150392 and GSE167028 were obtained from the Gene Expression Omnibus (GEO), including cardiomyocytes derived from human induced pluripotent stem cells infected with SARS-CoV-2 in vitro and GSE150392 from patients with myocarditis infected with SARS-CoV-2 and the GSE167028 gene expression dataset. Differentially expressed genes (DEGs) (adjusted P-Value <0.01 and |Log2 Fold Change| ≥2) in GSE150392 were assessed by NetworkAnalyst 3.0. Meanwhile, significant modular genes in GSE167028 were identified by weighted gene correlation network analysis (WGCNA) and overlapped with DEGs to obtain common genes. Functional enrichment analyses were performed by using the "clusterProfiler" package in the R software, and protein-protein interaction (PPI) networks were constructed on the STRING website (https://cn.string-db.org/). Critical genes were identified by the CytoHubba plugin of Cytoscape by 5 algorithms. Transcription factor-gene (TF-gene) and Transcription factor-microRibonucleic acid (TF-miRNA) coregulatory networks construction were performed by NetworkAnalyst 3.0 and displayed in Cytoscape. Finally, Drug Signatures Database (DSigDB) was used to probe drugs associated with COVID-19 Myocarditis. RESULTS: Totally 850 DEGs (including 449 up-regulated and 401 down-regulated genes) and 159 significant genes in turquoise modules were identified from GSE150392 and GSE167028, respectively. Functional enrichment analysis indicated that common genes were mainly enriched in biological processes such as cell cycle and ubiquitin-protein hydrolysis. 6 genes (CDK1, KIF20A, PBK, KIF2C, CDC20, UBE2C) were identified as critical genes. TF-gene interactions and TF-miRNA coregulatory network were constructed successfully. A total of 10 drugs, (such as Etoposide, Methotrexate, Troglitazone, etc) were considered as target drugs for COVID-19 Myocarditis. CONCLUSIONS: Through bioinformatics method analysis, this study provides a new perspective to explore the pathogenesis, gene regulatory networks and provide drug compounds as a reference for COVID-19 Myocarditis. It is worth highlighting that critical genes (CDK1, KIF20A, PBK, KIF2C, CDC20, UBE2C) may be potential biomarkers and treatment targets of COVID-19 Myocarditis for future study.

Mechanism of COVID-19-Related Proteins in Spinal Tuberculosis: Immune Dysregulation.

Purpose: The purpose of this article was to investigate the mechanism of immune dysregulation of COVID-19-related proteins in spinal tuberculosis (STB). Methods: Clinical data were collected to construct a nomogram model. C-index, calibration curve, ROC curve, and DCA curve were used to assess the predictive ability and accuracy of the model. Additionally, 10 intervertebral disc samples were collected for protein identification. Bioinformatics was used to analyze differentially expressed proteins (DEPs), including immune cells analysis, Gene Ontology (GO) and KEGG pathway enrichment analysis, and protein-protein interaction networks (PPI). Results: The nomogram predicted risk of STB ranging from 0.01 to 0.994. The C-index and AUC in the training set were 0.872 and 0.862, respectively. The results in the external validation set were consistent with the training set. Immune cells scores indicated that B cells naive in STB tissues were significantly lower than non-TB spinal tissues. Hub proteins were calculated by Degree, Closeness, and MCC methods. The main KEGG pathway included Coronavirus disease-COVID-19. There were 9 key proteins in the intersection of COVID-19-related proteins and hub proteins. There was a negative correlation between B cells naive and RPL19. COVID-19-related proteins were associated with immune genes. Conclusion: Lymphocytes were predictive factors for the diagnosis of STB. Immune cells showed low expression in STB. Nine COVID-19-related proteins were involved in STB mechanisms. These nine key proteins may suppress the immune mechanism of STB by regulating the expression of immune genes.

Immune-Related Protein Interaction Network in Severe COVID-19 Patients toward the Identification of Key Proteins and Drug Repurposing.

Coronavirus disease 2019 (COVID-19) is still an active global public health issue. Although vaccines and therapeutic options are available, some patients experience severe conditions and need critical care support. Hence, identifying key genes or proteins involved in immune-related severe COVID-19 is necessary to find or develop the targeted therapies. This study proposed a novel construction of an immune-related protein interaction network (IPIN) in severe cases with the use of a network diffusion technique on a human interactome network and transcriptomic data. Enrichment analysis revealed that the IPIN was mainly associated with antiviral, innate immune, apoptosis, cell division, and cell cycle regulation signaling pathways. Twenty-three proteins were identified as key proteins to find associated drugs. Finally, poly (I:C), mitomycin C, decitabine, gemcitabine, hydroxyurea, tamoxifen, and curcumin were the potential drugs interacting with the key proteins to heal severe COVID-19. In conclusion, IPIN can be a good representative network for the immune system that integrates the protein interaction network and transcriptomic data. Thus, the key proteins and target drugs in IPIN help to find a new treatment with the use of existing drugs to treat the disease apart from vaccination and conventional antiviral therapy.

Insights into potential mechanisms of asthma patients with COVID-19: A study based on the gene expression profiling of bronchoalveolar lavage fluid.

BACKGROUND: The 2019 novel coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a major challenge threatening the global healthcare system. Respiratory virus infection is the most common cause of asthma attacks, and thus COVID-19 may contribute to an increase in asthma exacerbations. However, the mechanisms of COVID-19/asthma comorbidity remain unclear. METHODS: The "Limma" package or "DESeq2" package was used to screen differentially expressed genes (DEGs). Alveolar lavage fluid datasets of COVID-19 and asthma were obtained from the GEO and GSV database. A series of analyses of common host factors for COVID-19 and asthma were conducted, including PPI network construction, module analysis, enrichment analysis, inference of the upstream pathway activity of host factors, tissue-specific analysis and drug candidate prediction. Finally, the key host factors were verified in the GSE152418 and GSE164805 datasets. RESULTS: 192 overlapping host factors were obtained by analyzing the intersection of asthma and COVID-19. FN1, UBA52, EEF1A1, ITGB1, XPO1, NPM1, EGR1, EIF4E, SRSF1, CCR5, PXN, IRF8 and DDX5 as host factors were tightly connected in the PPI network. Module analysis identified five modules with different biological functions and pathways. According to the degree values ranking in the PPI network, EEF1A1, EGR1, UBA52, DDX5 and IRF8 were considered as the key cohost factors for COVID-19 and asthma. The H2O2, VEGF, IL-1 and Wnt signaling pathways had the strongest activities in the upstream pathways. Tissue-specific enrichment analysis revealed the different expression levels of the five critical host factors. LY294002, wortmannin, PD98059 and heparin might have great potential to evolve into therapeutic drugs for COVID-19 and asthma comorbidity. Finally, the validation dataset confirmed that the expression of five key host factors were statistically significant among COVID-19 groups with different severity and healthy control subjects. CONCLUSIONS: This study constructed a network of common host factors between asthma and COVID-19 and predicted several drugs with therapeutic potential. Therefore, this study is likely to provide a reference for the management and treatment for COVID-19/asthma comorbidity.